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ABSTRACT 

We present the discovery of five new unbound hypervelocity stars (HVSs) in the outer Milky Way 
halo. Using a conservative estimate of Galactic escape velocity, our targeted spectroscopic survey has 
now identified 16 unbound HVSs as well as a comparable number of HVSs ejected on bound trajec- 
tories. A Galactic center origin for the HVSs is supported by their unbound velocities, the observed 
number of unbound stars, their stellar nature, their ejection time distribution, and their Galactic 
latitude and longitude distribution. Other proposed origins for the unbound HVSs, such as runaway 
ejections from the disk or dwarf galaxy tidal debris, cannot be reconciled with the observations. An 
intriguing result is the spatial anisotropy of HVSs on the sky, which possibly reflects an anisotropic 
potential in the central 10-100 pc region of the Galaxy. Further progress requires measurement of 
the spatial distribution of HVSs over the southern sky. Our survey also identifies seven B supergiants 
associated with known star-forming galaxies; the absence of B supergiants elsewhere in the survey 
implies there are no new star-forming galaxies in our survey footprint to a depth of 1-2 Mpc. 
Subject headings: Galaxy: halo — Galaxy: center — Galaxy: kinematics and dynamics — stars: 
early-type — galaxies: individual (M31, Sextans B) 



1. INTRODUCTION 

Unbound radial velocities distinguish HV Ss from other 
main sequence stars in the Galaxy. In iBrown et al.l 
(|2005[) we reported the first HVS: a short-lived 3 Mq 
star traveling with a Galactic rest frame velocity of 
700 ± 12 km s"'^, twice the Milky Way's escape velocity 
at the star's distance of ~100 kpc. The observed mo- 
tion of the star is comparable to the escape velocity from 
the surface of the star and is thus difficult to explain 
with stellar dynamics. The maximum ejectio n veloc- 
ity from stellar binary disruption mechanisms (jBlaauwl 
119611: iPoveda et all I1967D is limited to ~300 km s~\ 
for 3 Mf7^ stars (iLeonard &: Duncan 1988, 1990; Leonard 
1991lll993HTauris fc Takens 1998: Portegies Zwart 2000: 



Davies et all [20021: iGualandris et al.|[2005l ). 

There is overwhelming evidence for a 4 x 1 0^ Mp, mas- 
sive black hole (MBH) in the Galactic center ()Ghez et al.l 
I2008t iGillessen et"all I2009by The MBH sits in a vast 
crowd of stars, including sho rt-lived B stars with or bital 
periods as short as 15 years (jGillessen et al.l[2"009al ) and 
pericenter velocities as high as 12,000 km s~^ (=4% of 
the speed of light) (|Ghez et al.l 120050 . Three-body in- 
teractions between stars and the MBH are inevitable 
in this environment and will naturally r esult in th e 
ejection of unbound "hypervelocity stars" (jHillsl IT988I ) . 
This ejection process also works for a binary MBH 
(lYu fc Tremaind[200l . 

A MBH in the Galactic center must produce HVSs, and 
known HVSs fit the MBH ejection picture. The observed 
unbound velocities, stellar nature, ejection time distri- 
bution, and Galactic latitude and longitude distribution 
of HVSs all support a Galactic center origin. Known 
HVSs have the spectral types of B stars, the same spec- 
tral type as the stars observed orbiting the central MBH. 
In all cases where photometric variability (Fucntes et al.' 
2006f) or echelle spectroscopv (Przvbilla et al. 2008b. g: 



confirmed to be short-lived main sequence B stars. Even 
the observed number of HVSs is consiste nt with the the- 
oretically predi cted MBH ejection rate (|Yu fc Tremaind 
[2003: Pcrcts et al.l |2007| ) and with t he number of 
"S-stars" pres e ntly orbiting the MBH ([Peretsl [2009al : 
iBromlev et al.l I2012D . In the Hills three-body ex- 
change scenario, S-stars are th e former companions of 
HVSs (iGinsburg fc Loebl I2006D . a picture which may 
(jPerets et al.ll2009al ) ormav not ([Madigan et al.ll201l( ) 
be supported by the S-stars' eccentricity distribution. 
Proper motion measurements promise to more directly 
test the link be tween HVSs and the MBH in the Galac- 
tic center (e.g.. IBrown et al.]|2010b[ ). 

Here, we present the results of our on-going spectro- 
scopic survey on the 6.5m MMT telescope to find new 
HVSs. Our effective survey strategy is to target stars 
with the colors of ~3 Mq stars - stars bluer than halo 
blue horizontal branch (BHB) stars but redder than fore- 
ground white dwarfs - that should not exist at faint mag- 
nitudes unless they were ejected into the outer halo as 
HVSs. In previous papers we reported the discovery of 
14 unbound HVSs and a comparable population o f pos- 
sibly boun d HV S s over a fifth of the sky (Brown et al.l 
[2b06a b, 2Q07BI3, l2009al |O). Here we report 5 new un- 
bound HVSs; the distribution of HVS angular positions 
on the sky is significantly anisotropic. 

In §2 we describe our HVS survey strategy and spec- 
troscopic observations. In §3 we discuss the distribution 
of stars in the survey, and identify the bound and un- 
bound HVSs. In §4 we discuss the HVSs, and in §5 we 
discuss their spatial anisotropy. We conclude in §6. 



2. DATA 



2.1. 



Lopez-Morales fc Bonanosll2008() is available. HVSs are 



Target Selection Strategy 

We select HVS candidates with the magnitudes and 
colors of 2.5-4 Mq stars in the outer halo. 2.5-4 Mq stars 
are luminous and allow us to explore a large volume of 
space. Selecting these stars maximizes our contrast with 
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Fig. 1. — Color-color diagram showing the target seleetion for 
the unified HVS Survey (solid cyan line) and the observed can- 
didates (black squares). The Survey samples 17 < 3 < 20.25 to 
the left of the dashed line and 19 < g < 20.25 to the right of the 
dashed line as described in the text. Five new HVSs (solid ma- 
genta stars) and the fourteen previous HVSs (open magenta stars) 
scatter around the iGirardi et aU (2004) stellar evolution tracks for 
2.5-4 Mq main sequence stars (solid blue lines); numbers along 
the tracks indicate stellar mass. Average SDSS photometric un- 
certainties for a go = 19-5 star are indicated by the errorbars on 
the lower left. For reference, we plot the underlying distribution of 
stars in SDSS (black dots) and label the locus of white dwarfs and 
A-type/BHB stars. 



respect to the Milky Way's stellar population, because 
2.5-4 Mq stars are bluer than halo BHB stars and redder 
than disk white dwarfs. Targeting faint stars maximizes 
our efficiency for detecting HVSs because the density of 
normal Milky Way stars is very low in the outer halo. 
Finally, 2.5-4 Mq stars have main sequence lifetimes of 
a few X 10^ yr and thus should not exist at large distances 
unless they were ejected there. 

We select candidates based on de-reddened point 
spread function photometry from the Sloan Digital Sky 
Survey (SDSS, Aihara ct al. 2011). Each SDSS data re- 
lease changes its photometric calibration, however, which 
means that each data release contains candidates unique 
to that data release plus candidates in common with 
other data releases. Because we observe candidates taken 
from each subsequent data release (DR), our HVS Sur- 
vey is a concatenation of candidates selected from SDSS 
DR6, DR7, and DR8. 

We begin with a cut on reddening E{B — V) < 0.1 to 
ensure accurate colors. We exclude the small region of 
the SDSS between b < -1/5 + 50° and b>l/5- 50° to 
avoid excessive contamination from Galactic bulge stars. 
We also impose —0.5 < (r — i)o < to eliminate non- 
stellar objects such as quasars. 

Figure[T]illustrates our color selection. This color selec- 
tion unifies the orig inal HVS survey that t argets b righter 
and bluer objects (jBrown et all |2006a''b', '20 07bll3 ) and 
our more recent HVS survey that targets fainter and red- 



der objects (jBrown et al.ll2009al) . For reference, we plot 
in Figure [1] the underlying distribution of stars in SDSS 
and main sequence tracks for 2 M(^, 2.5 M^, 3 M p, and 
4 Mq solar metaUicity stars (jGirardi et all 120021 12004'). 
Our color selection region targets the sequence of 2.5 - 
4 Mq stars bounded by -0.4 < (g - r)o < (-0.43(m - 
g)o + 0.18) and (2.2(g - r)o + 1.1) < {u - g)o < 1.07. 

Our magnitude selection is 17 < go < 20.25 in the 
region {g — r)o < (— 0.43(u — g)o + 0.13). In the region 
(5 ~ '')o > (— 0.43(u — (7)0 + 0.13) we restrict ourselves to 
fainter magnitudes 18.6 -I- 10[(.g — r)o + 0.3] + [{u — g)o — 
(1.2(,g — r)o + 1-25)] < go < 20.25, or approximately 
19 <gQ < 20.25. 

Applying these selection criteria to the SDSS DR6, 
DR7, and DR8 photometric catalogs results in 874, 1097, 
and 1424 candidates, respectively. After removing ob- 
jects such as nearby galaxies and bright stars by visual 
inspection, 1509 unique candidates remain. We previ- 
ously published spectroscopic identifications and radial 
velocities for 609 of these candidates (plus a few hun- 
dred others that f all outside the presen t selection crite- 
ria) summarized in lBrown et al.l (|2010b[ ). SDSS provides 
spectroscopy for another 63 candidates. Thus there re- 
main 837 candidates to be observed; we report observa- 
tions of 497 of these candidates here. 

2.2. Spectroscopic Observations 

We obtained spectroscopy for 497 candidates at the 
6.5m MMT telescope during observing runs spanning 
December 2008 to October 2011. All observations 
were obtained with the Blue Channel spectrograph 
(| Schmidt et al.|[l989l ) using the 832 line mm~^ grating 
in second order and either a 1" or 1'.'25 slit. These set- 
tings provide wavelength coverage 3600 A to 4500 A and 
a spectral resolution of 1 - 1.2 A. All observations were 
obtained at the parallactic angle and were paired with 
comparison lamp exposures. 

Our goal was to obtain modest signal-to-noise (S/N) 
observations adequate for determining radial velocity. 
At g = 19 mag we typically used a 390 s exposure to 
obtain (S/N) ~ 7 pix~^ in the continuum. We pro- 
cessed the spectra in real-time to allow additional ob- 
servations of interesting candidates. We extracted the 
spectra using IRAFQ in the standard way and mea- 
sured rad ial velocities using th e cross-correlation package 
RVSAO (jKurtz fc Minklll998f) . The average statistical 
uncertainty of the radial velocity measurements is ±11 
km s~^. 

While we made an effort not to observe stars with ex- 
isting spectrosco py, SDSS re-o bserved 245 of our stars as 
part of SEGUE (JYannv et aLll2009D . We can use this in- 
dependent data set to verify our velocity calibration. The 
mean difference between our heliocentric radial velocities 
and elodiervfinal from SDSS DR8 is 2 ± 27 km s^^. The 
dispersion is consistent with the statistical uncertainties; 
2 km s^^ is our systematic uncertainty. 

Two epochs of spectroscopy also allow us to identify 
objects with variable radial velocity. Although no HVS 
exhibits a significant change in radial velocity, there are 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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eight probable white dwarfs that exhibit >4-cr changes 
in velocity between the two epochs. We expect that 
these objects are compact binaries, and we are pursu- 
ing follow-up observations to confirm their velocity vari- 
ability. To date, follow-up of low mass <0.25 Mq white 
dwarfs in the HVS Survey has resulted in the discovery of 
two dozen merging, double-degener ate binaries with or- 
bital p er iods as short as 12 minutes (iBrown et al.|[2010cl 
IMlallbl . [2QT2I: IKilic et al.l[20Tol l2011al lblld. l20'T2ir 

2.3. Survey Completeness and Spectroscopic 
Identifications 

With 1169 spectra in hand, the HVS survey is now 90% 
complete for DR6- and DR7-selected candidates and 43% 
complete for DR8-selected candidates. The remaining 
candidates are concentrated at faint magnitudes 19.75 < 
50 < 20.25 but are spread relatively evenly over the sky. 
This distribution of incompleteness is the result of many 
poor observing runs that prevented observations of the 
faintest objects in the survey. 

Of the 1169 objects: 955 (82%) are normal stars with 
late-B spectral types; 192 (16%) are white dwarfs; 15 
(1%) are quasars at z ~ 2.4; and 7 (1%) are B supergiant 
stars. The width and shape of the hydrogen Balmer lines 
provide a measure of surface gravity over our range of 
color (tempe r ature ). We use line indices described by 
IBrown et al.l (|2003f) to estimate spectral types, and we 
perform stellar atmosphere model fits for the late B-type 
obj ects using an upgrad e d vers ion of the code described 
by lAllende Prieto et al.l (|2006l ) . The white dwarfs are 
studied in a ser ies of pape r s on lo w mass white dwarfs 
that begin with IKilic et al.l (' 2007a, bV The quasars were 
reported previouslv (jBrown et al...2009al) . We discuss the 
B supergiants below, and focus the remainder of this 
paper on the 955 late B-type stars in the HVS survey. 

2.4. B Supergiants in the HVS Survey 

Our spectra reveal seven B supergiants in the ~10,000 
deg^ of sky sampled by our survey. All seven B super- 
giants are associated with known star-forming galaxies in 
the Local Group. The first two B supergiants belong to 
the Leo A dwarf. We used those supergiants plus others 
to make t he first stellar veloci ty dispersion measurement 
of Leo A (jBrownet al.ll2007aD . 

Four new B supergiants belong to M31 and 
trace its rotation curve. The stars are SDSS 

J003748.658-l-395402.57, SDSS J004108.023-l-401337.08, 
SDSS J004828. 734-1-423158.58, and SDSS 

J004855. 598-1-424629.99. The four stars have a 

mean apparent magnitude of 50 = 19.3 and thus are 
located at the distance of M31 ()Vilardell et al.l |2010( ) 



for a typical B9 lb luminosity of M„ 



-5.0. The stars 



have angular separations of 1° - 2° from the center of 
M31 and locations consistent with the M31 disk. The 
two northern stars have systemic velocities of -1-230 km 
s^^ and -1-250 km s^^ with respect to M31, and the 
two southern stars have systemic velocities of —126 
km s""'^ and —200 km s~^ with respect to M31. The 
velocities perfectly m atch the Hi rotation curve of M31 
(|Corbelli et al.|[2010D : the B supergiants clearly belong 
to the outer disk of M31. 

The final B supergiant, SDSS J095951. 180-1-052124.52, 
lies within 2f5 of the center of Sextans B. The star's 
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Fig. 2. — Minimum Galactic rest-frame velocity ii^j- distribution 
for the 955 late B-type stars in the HVS survey. The best-fit Gaus- 
sian (dashed line) has dispersion 107 km s~^, excluding the 39 stars 
with l^r/l > 275 km s"-"^. The significant absence of Vj.f < —275 
km s~^ stars demonstrates that the positive velocity outliers are 
short-lived. Stars with v^-f > 400 km s~^ are unbound. 

295 ± 4 km s~^ heliocentric radial velocity is identical 
with the velocity of Sextans B (Falco et al. 1999). Thus 
this B supergiant almost certainly belongs to Sextans 
B, and its existence is further evidence for on- going star 
formation in this Local Group dwar f. For a distance 
modulus of (m - M)o = 25.56 ± 0.10 (|Sakai et al.lfT997[ ) 
the go — 19.22 star must have Mg = —6.3, consistent 
with a B9 la supergiant. 

With absolute magnitudes of Mg — —5 to —6, the HVS 
Survey can detect B9 Ib/a supergiants to a depth of 1-2 
Mpc. The only known star-forming galaxy undetected in 
our survey footprint is Pegasus, a d warf that is known t o 
have an extremely weak blue plume (jMassey et al.1120071 ) . 
M33 also falls in our footprint and was detected by our 
color selection, however we excluded M33 stars from the 
observing list because of obvious extended emission sur- 
rounding the stars. The absence of B supergiants else- 
where in the HVS Survey implies that there are no new 
star-forming galaxies in our survey footprint to a depth 
of 1-2 Mpc. 

3. RESULTS 

We now discuss the properties of our 955 late B-type 
stars and identify the unbound HVSs. For reference. 
Gala ctic escape velocity at the solar circle is 500-600 km 
s-i (jSmithet al.|[2007f) and at 50 kpc it is 300-400 km 
s-i (IGnedin et al.ll2010l) . 

3.1. Radial Velocity Distribution 

Figure [2] plots the distribution of line-of-sight veloc- 
ities corrected to the Galactic rest-frame Vrf for the 
955 late B-type stars in the HVS Survey. We calcu- 
late rest frame velocities using the local standard of rest 
from .Schonrich et al., (201Ci) and a 250 km s~^ circular 
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rotat ion velocity (jReid et alj l2009t [McMillan fc BinnevI 
[2OT0h : 



D^/ 



w/ieHo + 1 1 ■ 1 COS / COS 6+262.24 sin ^ cos 6 + 7.25 sin 6. 

This definition of Vrf changes our previous estimates of 
Galactic rest-frame velocity by up to 30 km s~^. We note 
that Vrf represents the radial component of velocity in 
the Galactic rest frame, and is thus a lower limit to the 
true velocity of any HVS. 

The 915 survey stars with \vrf\ < 275 km s^^ have a 
± 4 km s~^ mean and a 107 ± 5 km s^^ dispersion, 
consistent with a stellar halo population. The velocity 
dis tribution, howeve r, is not exactly Gaussian (Figure 
Wl- IKing et al.l (120121) show that 10-17% of the HVS Sur- 
vey stars may belong to the Sagittarius stream, which 
explains the excess of stars in the range —50 to —150 km 
s~^ and at +100 km s~^. 

Stars with Vrf > +400 km s~^ are unbound and prob- 
able HVSs. We observe no stars moving towards us with 
Vrf < —400 km s~^, consistent with the picture that 
stars with Vrf > +400 km s~^ are ejected from the Milky 
Way. The most negative velocity star in o ur sample has 
Vrf = —359 ± 10 km s~^, consistent with iKenvon et all 
(j2008[) 's estimate of escape velocity at 50 kpc. If the un- 
bound stars are HVSs, then HVS ejection models show 
that we must also find a comparable numbe r of bound 
HVS s in our survey (Bromley et al. 2006: Kenvon et al.l 
12001 . 

Indeed, we observe a significant excess of stars around 
+300 km s^^ that are possibly bound HVSs. As an esti- 
mate of significance, there is less than a 10~^ probability 
of randomly drawing 18 stars with 275 < Vrf < 325 km 
s~^ from the tail of a Gaussian distribution with the ob- 
served parameters. Thus the excess of stars around +300 
km s~^ are probably not halo stars. Because the stars 
around +300 km s~^ are bound to the Milky Way, the 
absence of a comparable number of stars at —300 km 
s~^ demonstrates that the bound stars must have main- 
sequence lifeti mes less than their '- ^ l Gyr orbital turn- 
around times (iBrown et all l2007cl : iKollmeier fc GouJdl 
120071: lYufc Madaull2007[) . We simply do not see stars 
falling back onto the Galaxy with similar velocities. 
Given the observed colors and spectra, we conclude that 
the stars around +300 km s^^ are main-sequence ~3 Mq 
stars ejected into the outer halo. 

3.2. Unbound Stars 

We require distance estimates to determine whether a 
high velocity star is unbound. The observations demon- 
strate that high velocity stars are short-lived main se- 
quence stars. W e estimate lum inosities by comparing 
observed colors to IGirardi et al.l (,2002. .2004) stellar evo- 
lution tracks (see Figure [T]) for solar abundance main 
sequence stars. We use both {u — g)o and (g — r)o col- 
ors to estimate luminosity as described in IBrown et al.l 
(j2010b"). We estimate that the precision of our lumi- 
nosities is ±0.23 mag based on propagating color errors 
through the stellar evolution tracks with a Monte Carlo 
technique. We then calculate heliocentric distances using 
de-reddened apparent magnitudes, and convert to Galac- 
tocentric distances R assuming that the Sun is 8 kpc from 
the Galactic center. 
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Fig. 3. — Minimum rest-frame velocity vs. Galactocentric dis- 
tance R for all 37 stars with v^f > +275 km s"'^. Distances are 
estimated using IGirardi et al .' (12004') main sequence star tracks. 
Errorbars show the span of physicall y pos sible distance given 
by evolved BHB tracks l lDotter ct al] l200l) . Five new HVSs 
(solid magenta stars ) have velocities and distances exceeding the 
IGnedin et al.l II2010I ) escape velocity model (dashed line). Other 
possible HVSs (solid blue dots) have velocities and distances near 
the Kenyon ct al. (2008) escape velocity model (long dashed line). 
Previously identified HVSs are marked with open magenta stars. 
Isochrones of travel time from the Galactic center (dotted lines) 
are calculated assuming the observed minimum rest frame velocity 
Vrf is the full space motion of the stars. 

Systematics are the domin ant source of uncerta i nty in 
our distance estimates. The IGirardi et all ()2002l 120041) 
absolute magnitudes are 0.24 + 0.33 mag fainter than the 
SchalLer et al. (1992) absolute magnitudes that we used 
in earlier HVS survey papers. Thus our new luminosity 
estimates shift known HVSs to smaller distances. HVS7 
suffers the biggest change, going from 60 kpc to 37 kpc; 
echelle spectroscopy, however, establishe s that HVS 7 is 
a3/7 Mq Bp star with My ^ -1.05 ()Przvbilla et al.l 
l2008c[ ) and thus ~60 kpc distant. This example high- 
lights the uncertainty inherent in stellar evolution tracks. 
Stellar nature is another contributor to systematic uncer- 
tainty. We estimate a physical lower distance limit for 
our stars by assum i ng tha t they are evolved BHB stars. 
iDotter et all (|2007l . I2008D tracks show that hot, metal- 
poor BHB stars are 1.33 ± 0.69 mag fainter than a main 
sequence star of the same temperature in our survey. We 
indicate BHB distances with the horizontal error bars in 
Figure IH 

Figure O plots the resulting distribution of minimum 
Galactic rest frame velocity Vrf versus Galactocentric 
distance R for the 37 stars with Vrf > +275 km s~^. 
Identifying the unbound stars requires a potential model 
for the Galaxy. Because the Galactic potential is poorly 
constrained at large distances, we consider two recent 
models. 

IKenvon et al.l ()2008D construct a spherically symmetric 
potential to fit observed Milky Way mass measurements 
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TABLE 1 

HVS Survey Stars with v^f > +275 km s 



ID 



t>0 
(km s-i) 



'"rf_ 

(km s -*-) 



90 
(mag) 



(mag) 



Rgc 

(kpc) 



I 
(deg) 



b 
(deg) 



Catalog 



Rcf 



HVSs 



HVSl 


840 


673 


19.69 


+0.41 


76 


227.33 


+31.33 


SDSS J090744.99+024506.88 


1 


HVS2 


708 


718 


19.05 


+2.6 


26 


175.99 


+47.05 


US708 


2 


HVSS 


723 


520 


16.20 


-2.7 


62 


263.04 


-40.91 


HE0437-5439 


3 


HVS4 


611 


559 


18.31 


-0.66 


69 


194.76 


+42.56 


SDSS J091301. 01+305119.83 


4 


HVSS 


553 


664 


17.56 


+0.04 


36 


146.23 


+38.70 


SDSS J091759. 47+672238.35 


4 


HVS6 


626 


511 


18.97 


+0.10 


61 


243.12 


+59.56 


SDSS J110557.45+093439.47 


5 


HVS7 


529 


396 


17.64 


-0.19 


37 


263.83 


+57.95 


SDSS J113312. 12+010824.87 


5 


HVSS 


489 


393 


17.94 


+0.01 


46 


211.70 


+46.33 


SDSS J094214.03+200322.07 


6 


HVS9 


628 


461 


18.64 


+0.34 


47 


244.63 


+44.38 


SDSS J102137.08-005234.77 


6 


HVSIO 


478 


423 


19.22 


+0.21 


66 


249.93 


+75.72 


SDSS J120337.85+180250.35 


6 


HVS 12 


552 


408 


19.61 


+0.13 


82 


247.11 


+52.46 


SDSS J105009. 59+031550.67 


7 


HVS13 


575 


421 


20.02 


+0.09 


101 


251.65 


+50.64 


SDSS J105248. 30-000133.94 


7 


HVS14 


532 


396 


19.72 


-0.25 


96 


241.78 


+53.20 


SDSS J104401. 75+061139.02 


7 


HVS 17 


246 


442 


17.43 


-0.74 


50 


73.52 


+41.16 


SDSS J164156. 39+472346.12 




HVS18 


251 


473 


19.30 


+0.05 


75 


103.64 


-26.77 


SDSS J232904.94+330011.47 




HVS19 


597 


490 


20.06 


+0.35 


89 


256.05 


+63.74 


SDSS J113517.75+080201.49 




HVS20 


504 


381 


19.81 


+0.54 


69 


262.56 


+60.39 


SDSS J113637.13+033106.84 




HVS21 


355 


392 


19.73 


-0.23 


100 


165.26 


+56.11 


SDSS J103418. 25+481134.57 




Possible HVSs 


HVS 11 


477 


307 


19.58 


+0.68 


64 


238.77 


+40.63 


SDSS J095906. 47+000853.41 


7 




222 


292 


19.83 


+0.53 


77 


155.50 


+49.46 


SDSS J101359. 79+563111.65 


7 




496 


353 


18.86 


+0.28 


54 


256.27 


+54.55 


SDSS J111136. 44+005856.44 




HVS 15 


463 


322 


19.15 


-0.29 


78 


266.51 


+55.92 


SDSS J113341. 09-012114.25 


7 


HVS16 


434 


345 


19.33 


-0.03 


74 


285.86 


+67.38 


SDSS J122523. 40+052233.84 


7 










Possible Bound HVS 


5 








150 


321 


17.14 


+0.29 


27 


115.82 


-40.58 


SDSS J002810. 33+215809.66 






136 


311 


17.77 


+0.03 


40 


125.05 


-31.26 


SDSS J005956. 06+313439.29 






108 


276 


17.19 


+0.18 


30 


127.53 


-31.42 


SDSS J010948. 30+311727.66 






361 


288 


18.39 


-0.25 


61 


196.07 


+23.21 


SDSS J074950. 24+243841. 16 


7 




229 


297 


17.28 


+0.06 


34 


160.45 


+34.20 


SDSS J081828. 07+570922.07 


7 




306 


282 


18.08 


+0.14 


45 


186.30 


+42.16 


SDSS J090710. 07+365957.54 


7 




504 


339 


18.48 


+0.39 


44 


247.97 


+46.42 


SDSS J103357.26-011507.35 






448 


285 


19.23 


+0.14 


68 


250.71 


+47.87 


SDSS J104318. 29-013502. 51 






482 


309 


17.38 


+0.08 


30 


269.75 


+47.30 


SDSS J112255. 77-094734.92 






140 


291 


19.13 


+0.47 


58 


130.08 


+40.59 


SDSS J112359. 47+751807.73 






424 


290 


18.13 


-0.07 


44 


274.88 


+57.45 


SDSS J115245. 91-021116.21 


7 




228 


309 


17.49 


+0.06 


31 


65.34 


+72.37 


SDSS J140432. 38+352258.41 


7 




284 


289 


18.40 


+0.22 


40 


357.16 


+63.62 


SDSS J141723. 34+101245. 74 


7 




206 


278 


18.88 


-0.12 


58 


18.68 


+44.85 


SDSS J154806. 92+093423.93 






62 


292 


17.51 


+0.12 


29 


75.71 


+28.06 


SDSS J180050. 86+482424.63 






133 


290 


17.35 


+0.41 


25 


85.47 


-51.67 


SDSS J232229. 47+043651.45 





Referen ces. — ( 1) Brown ct al. (2005); (2) Hirsch ct al. (2005); (3) Edclmann et al.l l[2005h : (4) IBrown eFall 
(2006a); (5) IBrown et al. (2006b); (6) Brown ct al. (2007c); (7) Brown ct al. (2009;^) 



from 5 pc to 10^ pc. We estimate escape velocity from 
this model by dropping a test particle from rest at 250 
kpc and calculating its infall velocity. This calculation 
yields an escape velocity o f 360 km s~^ a t 50 k pc (long- 
dashed line in Figure [3). i Gnedin et al.l (|2010D measure 
the mass profile of the Milky Way using the velocity dis- 
persion of the (bound) HVS Survey stars. Assuming that 
the escape velocity is twice the circular velocity, we esti- 
mate an escape velocity of 400 km s~^ at 50 kpc (dashed 
line in Figure ^ . The escape velocity profiles differ in 
part because the IGnedin et al.l (20101) dark m atter halo 
is 1.6 X IQi^ M© while the lKenvon et al.l (|2008f ) dark mat- 
ter halo is 1 X 10^^ Mq. The difference in escape velocity 
profiles illustrates the present uncertainties. Because a 
larger halo mass is pr eferred by rec ent observations (e.g. , 
iPrzvbillaetaLJ |2Q10| ). we use the IGnedin et all (pOlO ) 
model to identify new unbound stars. 



We identify five new unbound HVSs wit h radial veloci- 
ties an d distances equaling or exceeding the lGnedin et al.l 
(|2010D escape velo city model, and one new possible 
HVS exceeding the iKenvon et al.l (|2008D escape veloc- 
ity model. Previous ly we identifie d 14 unbound HVSs 
on the basis of the IKenvon et al.l ([2008) model, how- 
ever our new Vrf and distance estimate s place HVSll 
slightly below the threshold of the Kenvon et al.l ()2008f ) 
model. Another tw o stars, HVS 15 and H VS16, are above 
the threshold of the lKenvon et al.l (I2008D ii iodel but well 
below the threshold of the IGnedin et al.l ()2010[) model. 
Therefore we re-classify these three objects as "possible 
HVSs." HVS7 and HVSS, on the other hand, are demon- 
strated :^4 M© main sequence B stars at R c^i^SO kpc 
(jLopez-Morales fc Bonanosll2008t iPrzvbilla et al.ll2008cD 
and thus are almost certainly HVSs. 

Table [1] summarizes the properties of all 37 HVS Sur- 
vey stars with Vrf > +275 km s^^. Columns include 
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HVS number, heliocentric radial velocity Vq, minimum 
Galactic rest-frame velocity Vrf (not a full space ve- 
locity), de-reddened SPSS q-band magnitude, absolute 
magnitude Mg from lGirardi et aD (|2002, 2004) main se- 
quence tracks, Galactocentric distance R, Galactic coor- 
dinates {l,b), and catalog identification. There are 16 
unbound HVSs, 5 possible HVSs, and 16 possibly bound 
HVSs. 

The velocity and distance estimates allow us to esti- 
mate ejection times. The dotted lines in Figure [3] in- 
dicate travel times from the Galactic center under the 
assumption that Vrf is the 3-dimensional velocity of our 
stars, a conservative assu mption even though the orbits 
should be nearly radial (jGnedin et al.|[2005l ). We find 
that HVSs are not clustered around a common ejection 
time, but have a range of ejection times spread continu- 
ously between 60 and 210 Myr. In other words, the ob- 
servations do not support a large burst of HVSs that one 
might expect for binary MBH in-spiral or dwarf galaxy 
tidal disruption models. Rather, the observations sup- 
port an on-going ejection process, consistent with the 
Hills three-body ejection model. 

4. FIVE NEW HYPERVELOCITY STARS 

We discuss the details of the five new unbound HVSs 
discovered in this portion of the survey. All of the HVSs 
except for HVS 17 have multiple epochs of observation. 
None of the HVSs exhibit significant radial velocity vari- 
ations, which indicates that the HVSs are unlikely to be 
compact binaries. 

SDSS J164156.391-h472346.12, hereafter HVS17, has 
a -f246 ± 9 km s~^ heliocentric radial velocity and a 
minimum velocity of Vrf = +443 in the Galactic frame. 
HVS17 is the hottest HVS in our survey with a B6 spec- 
tral type. HVS17 is also the brightest HVS in our sur- 
vey, g = 17.499 ± 0.015 mag, opening the possibility of 
follow-up echelle spectroscopy. A solar metallicity 4 Mq 
star with the temperatur e of HVS17 has Mg = —0.74 
(jGirardi et al.ll2002l [200l and a distance of i? = 50 kpc. 

SDSS J232904.947-h330011.47, hereafter HVS18, has 
a -1-251 ± 10 km s~^ heliocentric radial velocity and a 
minimum velocity of Vrf = -1-473 km s~^ in the Galac- 
tic frame. At b = -26.8°, HVS18 is the lowest latitude 
HVS in our survey. HVS 18 is also our first HVS discov- 
ery in the southern Galactic hem isphere. Unlike HVS3 , 
the 9 Mc:^ HVS near the LMC (jEdelmann et all 120051 : 
lPrzvbillaet"aIll2008b[ ). HVS18 is 17° from M31. Any 
association with M31 or its satellites is impossible. A so- 
lar metallicity 3 Mq star with the temperature of HVS18 
has Mg = 0.05, thus HVS18 has a distance of i? = 75 
kpc from the Milky Way and a distance of 700 kpc from 
M31. 

SDSS J113517.759-l-080201.49, hereafter HVS19, has 
a -1-597 ± 15 km s^^ heliocentric radial velocity and a 
minimum velocity of Vrf = -1-490 km s^^ in the Galactic 
frame. HVS19 is one of the cooler HVSs in our survey 
with an AO spectral type. A solar metallicity ~2.5 Mq 
star with the temperature of HVS19 has Mg — 0.35 
which locates it at i? = 90 kpc. Like many of the other 
HVSs, HVS19 is in the constellation of Leo. 

SDSS J113637.135-l-033106.84, hereafter HVS20, has a 
504 ± 12 km s~^ heliocentric radial velocity and a min- 
imum velocity of Vrf = -1-381 km s~^ in the Galactic 
frame. HVS20 sits on the escape velocity threshold of 



Gnedin et al.l (12010 
Kenvon et all (|200g 



but is well above the threshold of 
. HVS20 is similar to HVS19, only 
4?5 away; it is yet another HVS in the constellation Leo. 
HVS20 is also one of the coolest HVSs in our survey 
with an Al spectral type. With an estimated luminosity 
of Mg = 0.54, HVS20 sits at i? = 80 kpc and so shares 
a similar ejection time as HVS19. The other eight HVSs 
around Leo have ejection times that differ by up to 125 
Myr, however, implying that the HVSs in the direction 
of Leo do not all come from a single ejection event. 

SDSS J103418. 254+481134.57, hereafter HVS21, has 
a +355 ± 10 km s~^ heliocentric radial velocity and a 
minimum velocity of Vrf = +392 km s~^ in the Galactic 
frame. HVS21 has a B7 spectral type and, for a lu- 
minosity of Mg = —0.23, is located at R = 100 kpc. 
Like HVS20, HVS21 sits on escape velocity threshold of 
Gnedin et al.' (2010) but is well above the threshold of 
Kenvon et al. (2008). HVS21 is i n Ursa Major about 
11° from HVS2 ()Hirsch et al.ll2005[) . 

5. SPATIAL ANISOTROPY 

The spatial distribution of HVSs on the sky is intrigu- 
ing because it is almost certainly linked to their origin. In 
principle, the MB H in the Galactic cen ter can eject HVSs 
in all directions. iBrown et al.l (|2009bD demonstrate that 
unbound HVSs have a significantly anisotropic distribu- 
tion on the sky. Lower velocity bound HVSs have a more 
isotropic distribution. The new HVSs reported here sup- 
port these trends. 

Figure |4] plots the angular positions of all of the survey 
stars. The overall distribution of positions reflects the 
SDSS imaging footprint. HVSs are marked with solid 
magenta stars and possible HVSs are marked with solid 
blue dots. The stars have a similar distribution and so 
we consider the combined set of HVSs in an attempt to 
mitigate small number statistics. Because of the complex 
SDSS footprint and because of the HVS Survey incom- 
pleteness, we quantify the spatial anisotropy by compar- 
ing the distribution of HVSs against the distribution of 
observed stars from which they are drawn. Although this 
comparison is not perfect - the overdensity of stars arc- 
ing to the right of the North Galactic pole in Figure |4] 
is the Sgr stream - this approach provides a fair com- 
parison of positions and angular separations in the HVS 
survey. 

Figure [5] plots the cumulative Galactic longitude and 
latitude distributions of our survey stars and HVSs. 
Kolmogorov-Smirnov (K-S) tests find 0.013 and 0.415 
likelihoods that the combined set of HVSs are drawn 
from the same longitude and latitude distributions, re- 
spectively, as the survey stars. The longitude and lati- 
tude tests are not completely independent because of the 
complex boundaries of our survey region. The KS tests 
establish, however, that the HVS spatial anisotropy is 
primarily in Galactic longitude and not in Galactic lati- 
tude. 

Next we consider the distribution of angular separa- 
tions. Seven of the HVSs (33%) arc in the constellation 
Leo, thus we expect that the HVSs have an excess of 
small angular separations and an absence of large angu- 
lar separations compared to the survey stars. Calculat- 
ing the angular separations for all unique pairs of stars, 
a K-S test finds a 3 x 10"^* likelihood that the HVSs are 
drawn from the same distribution of angular separations 
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Fig. 4. — Polar projections, in Galactic coordinates, showing the 16 unbound HVSs (magenta stars), the 5 possible HVSs (blue circles) 
and 934 other stars (black dots) in the HVS survey over the southern (left) and northern (right) Galactic caps. 
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Fig. 5. — Cumulative distributions of Galactic I and b of the 
16 unbound HVSs (dashed magenta line), 21 HVSs and possible 
HVSs (dashed blue line), and the other stars (solid line) in the 
HVS survey. 

as the survey stars. We establish the significance of this 
hkehhood by taking 10,000 random draws of 21 unique 
stars from the survey: a hkehhood below 3 x 10~^^ occurs 
at random 0.3% of the time in our survey. The actual 
significance of the anisotropy depends on the method of 
calculation; our approach indicates a formal significance 
of 3-a. 

A number of models p ropose to explain the HVS spa- 
tial anisotropy. lAbadi et al. ( 20 0l) propo se a tidal debris 
origin for HVSs (see also iPiffi et al.ll2()lll ). However, this 
model is difficult to reconcile with the observations be- 



cause 1) known HVSs have ejection times that differ by 
up to 150 Myr, 2) no other unbo und tidal debris is ob- 
serve d in the same region of sky (jKoUmeier et al.l 120091 
l201Clf l. and 3) no dwarf galaxy in the Local Group has 
comparable velocity. An intermediate mass black hole 
in-spiral event in the Gala ctic center may produce a ring 
of HVSs around the sky (iGualandris et al.ll2005t iLevinI 
120061: iSesana et al.ll2006l:lBaumgardt et al.ll2006D but this 
model is also at variance with the observations. Theoret- 
ically predicted black hole in-spirals occur on timescales 
10-100 times shorter than the observed span of HVS ejec- 
tion times, thus multiple black hole in-spiral events with 
common alignment are necessary to account for the ob- 
served spatial distribution. 

Runaway stars, stars born in the disk and ejected by 
binary disruptions, arc another possible so urce of ap- 
parent HVSs (e.g. Gvaramadze et all l2009f ). Surveys 
of stars near the Galactic disk iden tify unbound stars 
ejected from the disk (Heber et al. 20^81: iPrzvbilla et al.l 
I2008at iTillich et al.) 120091: llrrgang et all I2010D as ex- 
pected in runaway star mechanisms (jBromlev et al.ll2009l : 
iGvaramadzd l2009h . Stars ejected from a rotating disk 
have a distinct spatial distribution: the fastest stars 
are those ejected in the direction of Galactic rotation, 
thus ~90% of >400 km s~ ^ runaways should be located 
|6| < 45° ()Bromlev et al.l [20091) . Moreover, >400 km 
s^^ runaways at i? = 50 kpc should be found over all 
Galactic longitudes. Thus the expected distribution of 
runaway latitudes and longitudes is inconsistent with the 
observed distribution of HVSs. 

The ejection rate of disk runaways is also inconsistent 
with the observed number of HVSs. The physical size 
of stars places a speed limit on stellar binary disruption 
mechanisms: a 3 Mq runaway with Vrf > -1-400 km s^^ 
at i? = 50 kpc requires the disruption of a compact bi- 
nary containing a massive star within its lifetime while 
avoiding a merger. An optimistic ejection rate for 3 Mq 
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runaways with >400 km s ^ at i? = 50 kpc is about 100 
to 100,000 times smalle r than the expected rate of HVSs 
electe d by the MBH (|Brown et al.ll2009al: iPerets fc Su"b^ 
|2012| ). Given that the HVSs are not clustered at low 
latitudes in our survey, we conclude that runaways are 
unlikely contaminants to the HVS sample. 

The anisotropic distribution of stars in the Galac- 
tic center provides another model for explaining the 
anisotropy of HVSs. If the central MBH is the origin of 
HVSs, and the stars orbiting near the MBH arc concen- 
tratcd in one or two disks (Levin & Bcloborodov 2003; 
ILuet al.l 120091: IBartko et ani2OO90 . then stars probably 
intera ct wi th the MBH in pr eferred planes. ILu et al.l 
(|2010D and lZhang et al.l (|2010[ ) demonstrate that the ob- 
served distribution of HVSs is consistent with two thin 
disk planes aligned with the present stellar disks in the 
Galactic center. Observed HVSs were ejected up to 200 
Myr ago, however. Although the Galactic center clearly 
contains anisotropic structure - including young stel- 
lar clusters like the Arches and Quintuplet clusters and 
gaseous disks, arms, and bars that may form young stars 
- there is no model that explains how Galactic center 
structure maintains a fixed orientation for 200 Myr. 

A final model is an anisotropic Galactic potential. 
HVSs are unique test particles because they start in the 
Galactic center and travel to infinity. The gravitational 
potential maps initial HVS ejection velocities to the ve- 
locities that we observe in the outer halo today. Because 
many HVSs are marginally unbound, a non-spherical po- 
tential can naturally explain why HVSs are found in pre- 
ferred directions on the sky: stars ejected along the long 
axis of the potential would have the fastest observed ve- 
locities; stars ejected along the minor axis of the poten- 
tial are decelerated more and so fewer will be observed 
as HVSs. Kenyon et al. (2008) show that the important 
quantity is the scale length of potential in the central 10- 
100 pc, the region where HVSs experience the most de- 
celeration. The anisotropic potential model can be tested 
with the all-sky distribution of HVSs. If there is rota- 
tion around the long axis of the potential, stars ejected 
along the long axis should be located in two clumps on 
opposite points in the sky. If there is rotation around the 
short axis, stars ejected along the long axis should ap- 
pear in ring of HVSs around the sky. There are currently 
no good constraints on the shape and orientation of the 
Galactic potential, but perhaps the HVSs will reveal it. 

6. CONCLUSIONS 

We report new results from our targeted HVS Survey, 
a spectroscopic survey of stars with ~3 Mq colors that 
should not exist at faint magnitudes unless they were 
ejected into the outer halo as HVSs. Our survey is 90% 
complete for SDSS DR6- and DR7-selected candidates 
and 43% complete for DR8-selected candidates. Spec- 
troscopy reveals that most of the candidates are normal 
late B-type stars in the halo; other objects include low 
mass white dwarf binaries and B supergiants in known 
star-forming galaxies. The absence of B supergiants else- 
where in the HVS Survey implies that there are no new 
star-forming galaxies to a depth of 1-2 Mpc within our 
survey footprint. 

The velocity distribution of the 955 late B-type stars in 
the survey is remarkable: 37 stars have Vrj > -1-275 km 
s~-^ and 3 stars have Vrf < —275 km s~^. The observa- 



tional signature of a HVS is its unbound velocity, which 
we determine by comparing observed radial velocities and 
distances to Galactic potential mo dels. There are 5 new 
unbound HVSs on the basis of the iGnedin et al.l ()2010l ) 
escape velocity model. The HVS Survey thus identifies 
16 unbound HVSs, 5 possible HVSs, and 16 possibly 
bound HVSs. 

A Galactic center origin for the HVSs is supported 
by their unbound velocities, the observed number of 
unbound stars, their stellar nature, their ejection time 
distribution, and their Galactic latitude and longitude 
distribution. Other proposed origins for the unbound 
HVSs, including runaway ejections from the disk or dwarf 
galaxy tidal debris, cannot easily be reconciled with the 
observations. Although not all unbound stars are neces- 
sarily HVSs, the MBH ejection origin provides the best 
explanation for the observed ~3 Mq unbound HVSs in 
our survey. 

HVSs are important because their properties are 
tied to the natu re and environment of the MBH 
that ej ects them (iPortegies Zwart et al.l l2006t iMerritt 
I2OOI IDemarque & Viranil 120071; 'Ginsburg & Loeb 
20071 : jGualandris & Portc gies_^Zwart 2007; Hansen 
2OO7I : ILu et al.l [2 007: Scs ana et al.l l2007allbl. 



2009^ lO'L earv fc Loeb. .2008tTPerets fc Alexandeii 
Percts 20 09allbl: 



2008; Lockmann 

2008, 2009; Chcn^laLJ 

2009 



Sherwin et al.l 20081 Svensson et al 
Baumgardti |200S ^ 'Lockmann et al 



20091 



Gualandris & Merritt 
2009; Sari et al. 



^. ._-, Homnan 2009; Madigan et al. , 

2010: Antonini et al. 2010 . 2011: Baruteau et al.l 120111 : 
Ginsburg fc PeretsI I2012D . HVS are also important 



probes of the da rk matter potenti al through which they 
move (IGnedin et al. 20 01 lYu fc M adau 2007: iWu et "all 
I200I iKenvon et al.ll2008nPeret7et al.ll2009bl) . 

An intriguing result of the HVS Survey is the spatial 
anisotropy of unbound HVSs on the sky. There is as yet 
no compelling explanation for the anisotropy. We spec- 
ulate that an anisotropic central potential may explain 
the observations. In the future, measuring the full space 
trajectori es of the HVSs prom ises to elucidate the origin 
of HVSs (Bro wn et al.ll2010al ) and better constrain the 
shape and orientation of the Galactic potential. Further 
progress requires measuring the spatial distribution of 
HVSs over the entir e sky. We look for ward to using the 
Skymapper Survey (jKeiler et al.ll2007[) to identify HVS 
candidates in coming years. 

We thank M. Alegria, A. Milone, and J. McAfee for 
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Observatory, a joint facility of the Smithsonian Institu- 
tion and the University of Arizona. This project makes 
use of data products from the Sloan Digital Sky Sur- 
vey, which is managed by the Astrophysical Research 
Consortium for the Participating Institutions. This re- 
search makes use of NASA's Astrophysics Data System 
Bibliographic Services. This work was supported by the 
Smithsonian Institution. 
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Note. — (This table is available in its entirety in machine-readable and Virtual Observa- 
tory forms in the online journal. A portion is shown here for guidance regarding its form and 
content.) 
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APPENDIX 

DATA TABLE 



Table O presents the sample of 955 late B-type stars in the HVS Survey studied here. For completeness, we provide 
both the 375 newly observed objects and the 580 previously reported objects. Columns include J2000 coordinates, 
SDSS g-hand apparent magnitude, our heliocentric velocity, velocity in the Galactic rest frame (Equation [1]), and 
Galactic longitude and latitude. 



